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I. INTRODUCTION 

CP violating phenomena receive strong and continu- 
ous interest because they provide a gateway to physics 
beyond the Standard Model (SM). While CP violation is 
observed in the quark sector and is currently believed 
to be predominantly stemming from the Kobayashi- 
Maskawa-mcchanism, further searches are being pursued 
to test this (SM-)picture where CP and flavor violation 
are intimately linked. No breakdown of CP symmetry 
has been seen so far among leptons, however, there has 
been great progress made in neutrino masses and mixing. 

We consider here lepton electric dipole moments 
(EDMs) £ E dm = dii-i^ia^ysF^^ [1], specifically 
for muons, as probes of CP and lepton flavor violation. 
Lepton EDMs in the SM appear first at four loop or- 
der, and are tiny, e.g., df M < 10 _38 ecm [2]. The current 
experimental limits are [3, 4] 

d e = (6.9 ±7.4) • 10~ 28 ecm, 

df, = (-0.1 ±0.9) • l(T 19 ecm. (1) 

While the bound on the electron EDM severely constrains 
flavor blind CP phases, lepton flavor violating couplings 
are subject to the constraints from the branching ratios 
of rare lepton decays. At 90 % C.L. [3, 5] 

B(n -> &i) < 1.2 • HT 11 , 

B(t -> ^7) < 4.4- 10~ 8 , 

B(t -> e 7 ) < 3.3 • 10~ 8 . (2) 

This situation can change significantly due to dedicated 
experimental initiatives in addition to the upcoming di- 
rect searches up to multi-TeV energies at the LHC: The 
MEG collaboration expects a reach in the fi — > 67 
branching ratio of order 10~ 13 in the next few years [6]; 
The bounds on the radiative tau decays can be improved 



at a possible future super flavor factory up to 2 • 10 -9 [7] 
(with 75ab _1 ); there is a recent proposal to improve the 
current limit on d M by three orders of magnitude, even 
as low as 5 • 10~ 25 ecm [8]. 

We ask here whether cZ M can be large, i.e., not many 
orders of magnitude below its current bound and if 
is large, what are the requirements and implications for 
flavor physics? 

We address these questions within the minimal su- 
persymmetric Standard Model (MSSM). New sources of 
lepton flavor violation arise from supersymmetry (susy) 
breaking contributions causing intergenerational slepton 
mixing, for earlier works, see e.g., [9-11], and [12] for 
flavor non-universal but diagonal effects. With current 
data, and depending on the mass scale of the susy spec- 
trum, for sizeable muon EDMs some fine-tuning is in- 
volved. We measure its amount within the general MSSM 
considering scenarios where CP violation is genuinely 
linked to lepton flavor violation only. In this frame- 
work, there is no model-independent connection between 
hadronic and leptonic CP violation, and we do not im- 
pose the CP constraints from the hadron sector. Specif- 
ically, nuclear effects in the extraction of the electron 
EDM from the Thallium EDM, see, e.g., [1, 13], and the 
2-loop effects from [13] are not included here. The direct 
link to weak CP violation highlights the importance of 
the muon EDM with respect to the nucleon ones. 

Models with gauge mediation being the dominant 
source of susy breaking but with an additional contribu- 
tion from Planck-scale gravity have recently been stud- 
ied for their non- minimal flavor properties [14, 15], but 
also the vacuum structure [16]. We work out the condi- 
tions for a large muon EDM for such a realistic hybrid 
model and argue that an observation in the range of the 
anticipated reach could be explained with very specific 
solutions to the flavor problem only. 

The plan of the paper is as follows: In Section II we 



introduce the basic MSSM parameters relevant for the 
evaluation of the leptonic EDMs and rare decays, and 
discuss the relation between these observables. A de- 
tailed numerical study in the general MSSM is presented 
in Section III, where we also access the fine-tuning re- 
quired for a large muon EDM. Models based on hybrid 
gauge-gravity mediation are investigated in Section IV. 
More conventions and formulae are given in the appendix. 

II. THE MUON EDM WITH FLAVOR 

In case of flavor blind CP violation only the muon 
EDM is constrained by the tight limit on the electron 
one, Eq. (I) to be below, at 90 % C.L., 

~ m^/m e d e < 3.9 • 10 _25 ecm, (3) 

where mi denotes the lepton masses. We are thus led to 
consider CP violation in flavor violation to obtain larger 
values of d^. For the purpose of this work we set all 
CP phases not related to flavor to zero. To ease the 
notation we use di interchangeably for both the EDM 
and its magnitude throughout this work. 

In Section II A we introduce the susy slepton flavor 
sector and define the requisite mass insertion parame- 
ters. Constraints from rare decays of lcptons Eq. (2) put 
constraints on the amount of flavor violation. In Section 
II B we discuss the interplay between the muon EDM 
and the rare decays. In Section II C we investigate the 
higgsino contributions. 

A. Slepton Flavor 

Genuine susy flavor violation enters through the struc- 
ture of the soft terms in generation space. This concerns 
the mass-squared matrix of the charged sleptons, which 
is given by 

M 2 ( M l L ML. \ (i) 

1 - \ Ml R ) ' (4) 

and which connects left-chiral II and right-chiral sleptons 
Ir as r N Mjj M lM, with the chiral projectors N, M — L, R. 
The 3x3 sub-matrices read as 

Mh tlJ = MIm + (F,D- terms), 

M LR,ij = VlA B ,ij - fJ,V 2 (Y e )ij, 

M RRtij = M% tij + {F,D- terms). (5) 

Here, i,j = 1,2,3 denote generational indices and Y e is 
the Yukawa matrix of the charged leptons. The Hig gs 
vacuum expectation values v±.2 obey v = y/v \ + uf ~ 
174 GeV and V2/V1 = tan/3. The fi parameter denotes 
the Higgs mass term from the MSSM superpotential. 

The sneutrino masses cause intergenerational mixing 
as well. In the presence of left-handed neutrinos only, 
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FIG. 1: Flavor and chirality flow of the leading diagram con- 
tributing to the muon EDM. A cross denotes a mass insertion. 
The photon is attached wherever possible. 

the mass-squared matrix is written as 

= M h 3 + {F,D- terms) . (6) 

In M RR LL and M~ we did not spell out explicitly the 
flavor diagonal F- and D-terms whose effect on the di- 
agonal entries is suppressed as mf /M 2 E and v 2 /Ml E , 
respectively. In the full numerical analysis these terms 
are included. 

To make contact with the low energy phenomenol- 
ogy, we need to evaluate the soft terms at the weak 
scale, m,z- Further, we go over from the flavor eigen- 
states to the basis where the leptons are mass eigenstates 
and the neutralino interactions are diagonal in genera- 
tion space. We denote the respective unitary matrix of 
slepton-type A — Ll,vl,E by Va, and the correspond- 
ing slepton soft terms by a tilde, M\ R = Vl l M^ r V e and 
M\ = V A M 2 A V\. From the latter the mass insertions can 
be read off as 

S&W = (M 2 A h/M 2 A , (7) 

where we introduced an average slepton mass Mj. The 
6 parameters Eq. (7) induce flavor changing neutral cur- 
rents (FCNCs) and if complex, CP violation. 

B. dfj, versus B(t — > fj,"f) 

To understand the interplay between a large muon 
EDM but small enough lepton flavor violating branch- 
ing ratios we employ the following approximations: We 
restrict ourselves to stau-smuon flavor mixing only, and 
neglect all lepton masses except for the one of the tau; 
we set to zero the flavor violating chirality-flipping cou- 
plings which are contained in the A^-terms of Eq. (5). 
Furthermore - note that this only matters for the rare de- 
cays - we assume that the leading contributions are due 
to photino exchange. We use the formulae from Ref. [17], 
which are obtained in the mass insertion approximation, 
i.e., for perturbative deltas sufficiently smaller than one. 

In this approximation, the leading contribution to d^ 
can be obtained at one-loop from 
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where y = M 2 /M|, and M; denotes the relevant gaugino 
mass, that is, here the photino. The loop function G\(y) 
obeys C?i(l) = 1/12, drops rapidly with increasing y, and 
is given in the appendix. 

To estimate the EDM from flavor we approximate 
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by its effective value 5^ 5^f E 5 E ^ , see also Fig. 1, 



where the bino contribution is shown. Then, for y — 1, 
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We allow for maximal CP phases of S 2 3 S E ^ . The fac- 
tor E ', see Eq. (7), parametrizes the left-right (LR) 
mixing of the staus, and is taken to be real here. In our 
analysis below we fix the value of S 3 ^ E , thereby linking 
the dependence on ^4 Bi33 ,/i,tan/3 and Af$. 

The r — > ^7 branching ratio can we written using the 
same approximations as 



B(r -> ^7) 



l^l 2 ) 



(10) 
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FIG. 2: The maximal value of the muon EDM from flavor as 
a function of the average slepton mass based on Eqs. (8)-(12), 
and y = 1. The three solid curves correspond to the current 
upper limit on B(t — > yy), Eq. (2), and are obtained for 
5^ E = 10~ 3 , 1CT 4 and 10" 2 (from top to bottom at M A = 
1000 GeV). The dashed line refers to the hypothetical limit 
B(t ->• yj) < 2 ■ 10~ 9 for 8^ E = 10~ 3 . The thick curves 
have mass insertions \5\ < 1. 
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We approximated and by their effective values 
633 and S E 2 S^ E , respectively. The loop function 
Gs(y) satisfies 63(1) = 1/40 and is given in the appendix. 

The maximal value for the muon EDM allowed by the 
upper limit on the branching ratio B(t — > /i7) max is then 
determined by [9] 



maxl^Sfal = B(t -s- ^7) max /( 2 «0- 



(12) 



The resulting reach is shown in Figs. 2 and 3. The 
curves in these figures obtained from extrapolations be- 
yond the validity of the mass insertion approximation 
(thin lines) arc expected to illustrate the qualitative fea- 
tures only. We used B(t — > [ivv) = 17.36% [3]. 

We learn the following: 

• Values of d^ as large as 0(lO~ 22 )ecm are consis- 
tent with current FCNC constraints. For this to 
happen it requires 0(1) intergenerational mixing. 
(The mass insertion approximation breaks down). 

• The maximal allowed value of the muon EDM given 
by Eq. (12) grows with increasing and y. 

• The muon EDM vanishes for vanishing stau LR- 
mixing 633 E in our approximation. Large LR- 
mixing, however, suppresses because of the en- 
hancement of the coefficient n. For our parameters 
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FIG. 3: The same as in Fig. 2 but with the sign of 5 33 L 
flipped, i.e., being negative. 



The sign of the stau chirality mixing matters: Neg- 
ative values of 6 33 L allow for larger EDMs by sup- 
pressing B(t — > fij), see Eq. (11), at the price of 
increased tuning. This is illustrated in Fig. 3. 



Within the mass insertion approximation, \5\ < 1, 
values of d^ up to 0(1O~ 23 ) ecm are possible if the 
sleptons have masses below a few hundred GeV, 



J 33 



10-( 2 ~ 3 ) and y < 0(1). 



The maximal value for d^ is proportional to the 
upper limit on B(t — > ^7). The anticipated future 
bound on B(t — > /ij) from super flavor factories of 
2 ■ 10~ 9 [7] has a significant impact on the maximal 
value of d,,. 
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FIG. 4: The maximal value of the muon EDM from flavor 
as a function of the average slepton mass based on [9] for 
tan ,3 = 3, Mi = M A and (i = 2 TeV. The three solid curves 
correspond to the current upper limit on B(r — > fry), Eq. (2), 
and are obtained for 5^ E = 10~ 2 ,10" 3 and 10" 4 (from top 
to bottom at M A = 1000 GeV). The dashed line refers to the 
hypothetical limit B{r -> hj) < 2 ■ 10~ 9 for 8^ E = 10 -3 . 
The thick curves have mass insertions |<5| < 1. 
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FIG. 5: The same as in Fig. 4 but with /i = M A , hence lighter 
and more important higgsinos. 



C. Including Higgsinos 

Since we neglect the mass of the muon, the higgsinos 
(and winos) do not contribute to to the order we are 
working, and Fig. 1 represents the leading contribution 
to the EDM. However, higgsino contributions have been 
found to be of importance in the calculation of B(t — > 
/ry) [9] . More specifically, for a light /i term /x ~ Mi , M 2 
the lightest neutralino has a substantial higgsino fraction 
and the photino-only approximation in the rare r decays 
Eq. (10) receives large corrections. Consequently, the 
bounds on flavor violation change in the presence of the 
higgsinos, which then affects the maximal value of 
from flavor. 

The higgsino contribution to the r — > /17 amplitude 
involves directly the t Yukawa coupling y T oc tan /3. 
Hence, keeping 8 a £ E fixed while enhancing tan/3, the 
higgsino loops get scaled up with respect to the bino 
ones and the maximal value of the muon EDM drops 



with respect to the analysis of the previous section. The 
same holds if the A-term dominates the stau LR-mixing, 
633 E ~ vAe 33/ tan/3. If the /i term is large and dom- 
inates, S^ E ~ —\xm T tan /3, then all bino and higgsino 
contributions grow with tan /3. Since the EDM depends 
only linearly on 5^ , whereas B(t — > /17) does quadrat- 
ically, the maximal value of the EDM drops with increas- 
ing tan j3 also in this case. Therefore, large values of d^ 
favor a low tan /3. 

In Fig. 4 the maximal muon EDM as a function of the 
common slepton (smuon) mass is presented for the same 
values of <5 33 L as in Fig. 2, however with the higgsino 
and wino contributions in B(t — > jj/j) included. We use 
tan/3 = 3, /z = 2 TeV, M x = M A and M 2 = {gl/gl)Mx. 
As expected, a reasonably high value of \i suppresses the 
contribution from the higgsino loops, and, in general, the 
increase in B{r — > /X7) is small. Therefore, the photino 
analysis of Sec. II B gives a good approximation of the 
maximum value of the muon EDM. However, for low val- 
ues of E the leading bino graphs in the r — > fJ/y am- 
plitude, which are proportional to 5^ , are suppressed, 
and external chirality flip graphs oc m T , which are usually 
sub-leading, become more important at the low slepton 
mass range. Therefore, a cancellation against the hig- 
gsino graphs occurs, which can be seen in Fig. 4 as a 
spike in the allowed maximal EDM curves. For larger 
values of S^ E the cancellation does not happen, since 
the external flip contribution can not compete with the 
leading contributions. 

In Fig. 5 we show the maximal muon EDM for an hig- 
gsino mass parameter /j, which equals the slepton mass 
M^. Consequently, the smaller /i term enhances the hig- 
gsino mediated loop contributions in B(t — > fj/y) and the 
maximal EDM is not as large as in Figs. 2 or 4. With 
increasing slepton mass also fi increases, and the charac- 
teristics of Figs. 4 and 5 become alike. 

Here, we assumed that 6%% is equal to 8^ ■ Typically 
the sneutrino loop gives the largest single contribution to 
the branching ratio, unless the higgsinos are decoupled, 
in which case the bino loops dominate. 

To summarize, the inclusion of higgsino contributions 
modifies the correlation between the maximal value of 
the muon EDM from flavor and B(t — > fij) given in the 
previous section for low values of the \x term or when 
cancellations occur. Values of d M at order 10 -22 ecm are 
possible with current B(r — \ /X7) constraints for order one 
flavor mixings. Values of order 10 -23 ecm can be reached 
within the mass insertion approximation, \b\ < 1. 



III. NUMERICAL ANALYSIS IN THE 
GENERAL MSSM 

We study the muon EDM in the general MSSM ne- 
glecting CP phases unrelated to lepton flavor mixing. 
Formulae used are given in the appendix. 
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A. Input and Constraints on Susy Parameters 



We generate sets of data points by randomly sampling 
susy parameters at the weak scale. We consider two sce- 
narios, a "light" and a "heavy" one with their respective 
parameter ranges defined in Table I. 





flavor diagonal 


flavor off-diagonal 


light 

heavy 


- 1 TeV 
3-5 TeV 


- f 00 GeV 
0- 3 TeV 



TABLE I: The sampling ranges for the flavor diagonal and off- 
diagonal mass parameters in the light and the heavy scenario. 



As can be seen from Table f, we use different ranges 
for the flavor diagonal and off-diagonal mass parameters. 
The former contain the gaugino masses Mi and M2, the /j, 
term, its susy breaking companion, B^, and the diagonal 



entries in the slepton soft terms * / M 2 



For simplicity we restrict our analysis to degenerate first 
and second generation masses n = M^,™ 22 . Di- 

agonal A-terms are assumed to follow the corresponding 
lepton Yukawa couplings Ae,u — A$yi. Furthermore, 
we include slepton flavor mixing between the second and 
third generation only. The respective off-diagonal param- 
eters are therefore A E ^3, ^£,32, 



M l,23 and 



-^1.23- 



We allow for an order of magnitude suppression of the 
off-diagonal soft terms with respect to the diagonal en- 
tries in the light scenario, whereas in the heavy scenario 
we allow the off-diagonal terms to be of the same order of 
magnitude as the diagonal ones. In both scenarios tan/3 
is sampled over the range 2 - 50. 

For each generated point we check the scalar potential 
against unboundedness from below, and the existence of 
charge or color breaking minima [19]. These bounds con- 
strain the off-diagonal soft contributions of the A-terms. 
The soft masses of the Higgs fields are solved for with the 
vacuum condition dV/dcj) = 0. The resulting spectrum 
at each generated point is checked against experimental 
constraints from direct searches [3] . To generate a muon 
EDM we introduce random CP phases in all flavor off- 
diagonal soft terms between and 2ir. Since we assume 
here the \x — e and r — e slepton mixing terms and all 
flavor blind CP violation to be zero, the experimental 
constraint from the electron EDM is automatically ful- 
filled. Contributions to fi — > ej decays from sneutrino 
loops are also suppressed and neglected. 

For each sample point, we calculate the contribution to 
the muon anomalous magnetic moment, Aa M , the muon 
electric dipole moment d M , and the branching ratio for the 
decay r — > ^7. Aa M is constrained to be less than f CT 9 
corresponding to roughly two sigma of the experimental 
and theoretical uncertainty [18]. The bound on Aa^ is 
always fulfilled in the heavy scenario. 
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FIG. 6: The distribution of the muon EDM d M vs. the branch- 
ing ratio B(t — > ^17). Shown is the initial distribution in the 
light scenario (circles) and the heavy scenario (squares), as 
well as the distribution after the random walk process in the 
light (triangles up) and the heavy scenario (triangles down) . 



B. Random Walk Analysis 

The interplay between the muon EDM and the r — > fij 
branching ratio for our sampled points can be seen in 
Fig. 6. For both the light and heavy spectrum the value 
of the EDM for points satisfying B(t —> /17) < 10~ 8 is 
typically less than 10 -24 ecm (squares and circles). This 
value is, however, not a hard upper bound on the muon 
EDM, because it stems from the stochastic nature of our 
parameter selection. Indeed, and as we will see, one can 
find parameters such that the constraint on the branching 
ratio is fulfilled and is large. The question then is how 
fine-tuned such points are. 

To obtain points with > (10~ 22 — 1CP 24 ) ecm obey- 
ing the r — » /i7 constraint we use a random walk process. 
Starting from a point p in parameter space, we randomly 
vary the parameters by a small amount, p — > p' = p + 5p, 
and check for all the constraints as well as for a larger 
value of the muon EDM. ff p' passes the test, the pro- 
cedure is repeated with p' as the new starting point. If 
we start out with a point p that has B{t — > ^7) > ICT 8 , 
we additionally require a decrease in the branching ratio 
before accepting p' . The requirements on the branching 
ratio and the muon EDM thus drive a point into the de- 
sired direction. Unless we exceed a predetermined num- 
ber of failed attempts, this procedure eventually yields a 
point with B(t —> /17) < 1 0~ 8 and a large value of dp as 
requested, after which the random walk is terminated. 

In Fig. 6 we show the result of this random walk in 
the heavy and the light scenario. We see a tendency of 
points that start out with very low branching ratios to 
migrate to larger ones. This is expected, as the branch- 
ing ratio has similar contributions as the EDM, which 
is being forced to increasing values. It should be noted 
that for the fixed number of iterations, we find that in 
the heavy scenario only about 2% of the original points 
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this fraction is 35%. 



ecm, whereas in the light scenario 



(> 




FIG. 7: The distribution of tan /3 before (solid line) and after 
(dashed line) the random walk process in the light scenario 
(inclusive data set). 



For each scenario we consider in the following two 
data sets, an inclusive one containing all the points that 
went through the random walk routine and an exclusive 
one containing only points that actually pass the EDM 
constraint eL > 10~ 22 ccm for the light scenario and 
d[i > 0.5 • 10 -22 ecm for the heavy scenario. 

The random walk process affects the distribution of 
the parameters. E.g., in the light scenario there is a 
transition from a near flat distribution of tan j3 to one 
which is peaked at small values, see Fig. 7. The latter 
behavior is in agreement with the findings of Sec. II C. In 
the next section we use this selection effect to construct 
a measure of fine-tuning. 



C. Measures of Fine- Tuning 

In order to better understand the fine-tuning needed 
for achieving an experimentally interesting c? M while re- 
specting the experimental B(t — > /17) constraint, we first 

L / R 

note that the subamplitudes in a 2 g , see Eq. (A9), cor- 
responding to neutralino and chargino 1-loop Feynman 
diagrams, need to be small or cancel each other to great 
accuracy. We find that, however, there is usually only 
one subamplitude which contributes a dominant part of 
the branching ratio. When we increase the EDM via our 
random walk routine this disparity is lessened, but all 
individual subamplitudes remain small enough for there 
to be no need for large cancellations. 

Similarly, Eq. (A7) suggests that the imaginary part 

R / L 

Of 

°22 needs to be as large as possible for a large EDM. 
Since in our framework there is only an imaginary part 
in the flavor off-diagonal 2-3 mixing, the subamplitudes 
in the initial data set exhibit a dominance of the real 
part over the imaginary part for all s. A priori, one 
would expect this relation to be reversed by the require- 
ment to maximize d M but instead we end up with a sam- 
ple, where the real and imaginary parts are within one 
order of magnitude of each other in the light scenario, and 
only slightly more dispersed in the heavy scenario. This 



is sufficient for generating a large EDM since the overall 
sizes of the subamplitudes have increased. We conclude 
that measures of fine-tuning which rely on looking at can- 
cellations between different subamplitudes in magnitude 
or in imaginary part are not useful. 

Another way of measuring fine-tuning is to look at the 
local structure of parameter space and see how sensitive 
observables are to variations of the parameters [23]. We 
find that this measure of fine-tuning decreases for points 
when they go through our random walk process. This 
is expected since this method of measuring fine-tuning 
effectively looks at the slope in parameter space of a given 
observable, and in our case we drive these observables to 
their extremum values, i.e., to a region where the slope 
vanishes. 

To discuss fine-tuning we instead consider the width 
of a parameter's distribution and how the random walk 
process changes that width (see Fig. 7). From several 
such changes in the width we then construct a relative 
"volume" of parameter space that is preferred by points 
that go through our random walk process. We do this 
by comparing the number of fixed and equi-distant bins 
in a parameter X necessary to account for a given frac- 
tion fr of points before, iVbeforej and after, A a f tcr , the 
random walk process. This way, the bin number mea- 
sures the width of a (normalized) distribution. The ra- 
tio, r x = (iVafter/A^before), gives us a measure of how 
much the range of the parameter X has shrunk (peaked) , 
rx < 1, or expanded (flattened), rx > 1. Smaller 
(larger) volume fractions rx then correspond to larger 
(smaller) tuning. For the tan /3-distributions shown in 
Fig. 7 we obtain r tan/3 = 0.25 (for fr = 0.8). 

We consider the tuning of several parameters, specif- 
ically, the CP phases of the four flavor off-diagonal soft 
terms, masses and tan/3. 

Before proceeding, we note a couple of caveats. First, 
while random walking the range of some parameters can 
expand beyond the range to which we have initially con- 
strained them, e.g., the diagonal soft terms in the heavy 
scenario are sampled between 3 and 5 TeV, and the ran- 
dom walk process moves and spreads this range out. 
Since this is an artificial diffusion effect (and we can ef- 
fectively choose rx to be whatever we want by adjusting 
the initial sampling range) we exclude such parameters 
directly from the measure of fine-tuning. Secondly, even 
though a parameters distribution remains unchanged, it 
may become correlated with another parameter. An ex- 
ample for this are the complex phases of the off-diagonal 
soft mass terms. Let 2<fix 23 = arg(Afjf 23 ) for X = L,E. 
While wc find for each 4>l 23 and 4>e 23 separately the ratio 
r ~ 1, their difference A</> = 4>l 23 ~ 4>e 23 has r^ ~ 0.47 
with sharp peaks around 7r/4 and 37r/4 (heavy scenario), 
as can be seen in Fig. 8. This correlation can be under- 
stood from Eq. 9, with \^{8%£ 8%£ E $2$) maximized. 

We investigate the phase differences of M x 23 , X = 
L,E and Ae,23, ^£,32, and assign the product of the in- 



o - - — 

<t>L 23 - <t>E 23 

FIG. 8: The distribution of 0l 23 — ^£23 before (solid line) 
and after (dashed line) the random walk process in the heavy 
scenario (exclusive data set). 



dividual volume changes as 

d0 = IJrA</> ( . (13) 

i 

Furthermore, we consider various ratios of mass param- 
eters. We examine the following ratios containing off- 
diagonal entries: 

\Mj 23 \ \M%^\ \A E , 23 \ \A E &\ 

Ml 33 ' Miss ' A ' A ■ 1 j 

The product of their individual volume fractions is de- 
noted by d a d- We look at some general ratios of param- 
eters in the slepton sector and the mass parameters Mi, 
M 2 , fi: 




The product of their individual volume fractions is de- 
noted by da- We then define a total change of volume of 
parameter space, dy, as 

d v = d^ ■ d od ■ d G ■ r ta „/3. (16) 

One would expect that the volume changes dx depend 
on the value of the fraction of points, fr. In Fig. 9 we 
show this dependence for dv ■ As can be seen, dy is very 
stable over a large range of fr. The unstable behavior 
for values of fr > 0.8 is due to the way statistical out- 
liers affect our analysis. A more sophisticated approach 
to measuring the width of parameter distributions by fit- 
ting, e.g., Gaussian curves could improve the stability. 
We also verify that the values of the rj's (and conse- 
quently dx's) remain stable under variations of spurious 
quantities such as changes in data binning. 



oa 0.2 0.4 0.6 0.8 if> 

fr 

FIG. 9: The parameter space volume fraction dv Eq. (16) 
as a function of fr in the light inclusive (solid line), light ex- 
clusive (dotted line), heavy inclusive (dashed line), and heavy 
exclusive (dash-dotted line) scenario. 



D. Discussion 

The various volume fractions for the two scenarios and 
data sets are shown in Table II for fr = 0.8. Comparing 
the different entries, we get an idea about fine-tuning 
caused by the random walk process. 

Notably, the effect of the CP phases becoming corre- 
lated (as shown in Fig. 8) emerges in the exclusive heavy 
scenario only. Also we see that a tuning of tan /3 occurs 
in the light scenario only (see Fig. 7). 





Inclusive 


Exclusive 




Light Heavy 


Light Heavy 


?"tan/3 


0.25 0.96 


0.16 0.95 


dod 


1.3 0.88 


2.9 10 


dj, 


0.98 0.96 


0.98 0.49 


dG 


1.8 69 


0.39 33 


dv 


2.2 59 


1.1 170 



TABLE II: The volume fractions of parameter space after the 
random walk for fr = 0.8 for the light and the heavy scenario. 
Smaller numbers correspond to larger tuning. 

The ratios of the off-diagonal to diagonal mass ele- 
ments which we have collected into d d, have values close 
to unity in the inclusive set but spread out, by up to an 
order of magnitude in the case of heavy scenario, when 
we look at the exclusive set. This means that the more 
fine-tuned the points get, the wider the range of accept- 
able mass ratios is. These ratios enter the loop contribu- 
tions of the branching ratio and EDM and their increase 
in size, which leads to the spreading out of the distribu- 
tion, is in line with the observations from the beginning 
of Sec. IIIC. 

Finally, we see that for the set of flavor diagonal mass 
ratios in dc, there is a fine-tuning effect in the light sce- 
nario that emerges for the exclusive set, whereas in the 
heavy scenario there is instead again a spreading out. 

It is apparent that the ranges for the mass ratios are 
more constrained in the light scenario. In the heavy see- 



nario, the individual loop contributions (which partly de- 
pend on these mass ratios) start out much smaller than 
in the light scenario and thus there is more "room" to 
spread out. The phases (f>L 23 and 4>e 23 becoming corre- 
lated is also a sign of insufficient size of the loop contri- 
butions. 

It needs to be emphasized that the change in parameter 
space volume is a relative measure of the amount of fine- 
tuning and useful for comparisons between different data 
sets. It is not an absolute measure of the ease of fine- 
tuning. Even comparing different rows in Table II is not 
trivial, e.g., fine-tuning the phases can have a greater or 
lesser effect on the EDM than fine-tuning the mass ratios. 

We find that both heavy and light scenarios appear 
to be not particularly tuned to give a large value of the 
muon EDM while respecting existing constraints. While 
most of the volume fractions dx in the light scenario are 
smaller than in the heavy one, from looking at Fig. 6 and 
noting that a large percentage (~98%) of points in the 
heavy scenario can not be tuned to high enough values of 
dfj,, we conclude that the heavy scenario is quantitatively 
harder to fine-tune. 



IV. A FLAVORED EDM IN HYBRID MODELS 

After analyzing the muon EDM within the general 
MSSM, we now consider the situation in an explicit 
model with a hybrid mechanism of susy breaking, gauge- 
gravity mediation. These models have recently received 
attention because they provide insights into the nature 
of flavor breaking due to the appearance of flavor struc- 
tures beyond the Yukawa couplings [14, 15, 21, 22]. Vi- 
able gauge-gravity models have the dominant source of 
susy breaking from gauge mediation, which is flavor-blind 
at the scale of mediation. Additional contributions arise 
from Planck-scale gravity, which generically affect flavor 
physics. We assume here that the latter is controlled by 
Froggatt-Nielsen (FN) symmetries [24], generating also 
the Yukawa matrices. 

At the scale of gauge mediation, mjj, we parametrize 
the soft terms of the sleptons as [14] 

Ml(m M ) = m 2 L (l + x 2 X L ), 

M 2 E (m M ) = m%(l + Xl X R ). (17) 

The new flavor structure from gravity with respect to 
gauge mediation is encoded in the hermitean matrices 
Xl,r, following from the charge assignments of the FN- 
symmetry. The factors x\ and x 2 implicitly defined by 
Eq. (17) quantify the relative size of the gravity versus 
the gauge contribution to the SU{2) singlet and doublet 
soft masses, respectively. In this way, the Xi are mea- 
sures of the separation between the messenger and the 
Planck scale, mp\. In minimal gauge mediation with Nm 
messengers (see, e.g., [26]), 

/ m M \ 2 ( 4tt \ 2 Cl 5 2 

[—, c t?— , ci = -,c 2 = -, (18) 

\mpi/ \ai{m M ) J % 6 3 



where we assumed that the F-terms from mpi couple also 
to gauge mediation. The ratio x\jx 2 as in Eq. (18) is of 
order one, and increases for lower values of the messenger 
scale, see Fig. 10. 




1Q 6 1Q 8 1() 10 1Q 12 1Q 14 10 16 

m M [GeV] 



FIG. 10: The ratio x\/x2 in minimal gauge mediation Eq. (18) 
as a function of the messenger scale. 

In Section IV A we give formulae for the flavor violating 
low energy parameters in gauge-gravity models. A flavor 
model with explicit FN charges is studied in Section IV B. 
We investigate the impact of trilinear A-terms in Section 
IVC. 



A. Hybrid Gauge-Gravity 

Starting from the soft slepton masses at the scale rriM 
of the form given in Eq. (17) and solving the MSSM 
renormalization group (RG) equations [25] , the soft slep- 
ton masses at the weak scale can be written in the fol- 
lowing approximate form: 

M 2 L {m z ) ~ m£(r 2 l + c L Y e Y^ + x 2 X L ), 
M 2 (m z ) ~ mitnl + c^Ye + XiXfl). (19) 

In the spirit of the FN flavor models, Eq. (19) and the 
following estimations are only accurate up to numbers of 
order one, which we indicate by " ~" . The running of the 
gravity- induced terms is such an order one effect. The 
coefficients cl,ce are negative, obey \cl\ < \ce\ and are 
of the order ~ l/(167r 2 ) lnfrnjiz/m^), which are at most 
of order one for tum near the GUT scale. 

The Ti coefficients include the leading RG correction 
to the flavor diagonal elements. Neglecting contributions 
from the sub-leading gauge couplings, one obtains an an- 
alytical expression for j-j at one-loop order (sec, e.g., [26]): 

r, = l + ^(r"\4^)^. (20) 
CjTT \Jin( mz ) af{m M ) J mf(m M ) 

Here, Mi denotes the respective gaugino mass and m 2 
equals rh 2 L for i = 2 and m 2 E for i = 1. In messenger mod- 
els of gauge mediation, the ratio M 2 /m 2 is determined 
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by a simple formula at the scale of mediation including 
only the leading gauge contributions: 



Mf(m M ) 
mf(m M ) 



= CiN 



M ■ 



(21) 



The RG parameters r, from Eqs. (20) and (21) obey r 2 > 
r*i > 1 and are of order one for not too many messengers. 
Since r.- L > 1, the RG effect strengthens the flavor blind 
entries hence reduces the flavor violation. The properties 
of the ri can be different in general gauge mediation [27] 
which does not exclude negative soft sfermion masses- 
squared at the mediation scale. 

In the basis where the leptons are mass eigenstates 
and the neutralino interactions are diagonal in generation 
space we obtain: 



Ml L (m z ) 
M 2 VL {m z ) 
M 2 E {m z ) 



m\{r 2 \ + c L D 2 e + x 2 V Ll X l vI l ), 
m 2 L (r 2 l + c L C/ t D 2 e U + x 2 V Vh X L Vj L ) , 
m|(nl + c E D 2 e + x x V E X R vl), (22) 



where U = Vl l V$ l denotes the Maki-Nakagawa-Sakata 

(MNS) matrix and D e = diag(y e , y M , y T ) = V Ll Y £ V e . 
Note that the parametrization Eq. (22) holds beyond 
minimal gauge mediation with perturbative messengers. 

The mass insertions of the charged sleptons are then 
given as (i ^ j) 



X2 
> 2 



(V LL X L V! 



■ILL 
n 



(V E X R VI) IV (23) 



Here, we neglected contributions in the denominators 
from the tau-yukawa of the order y 2 ~ 10~ 4 tan 2 j3 with 
respect to the r, for mixings involving the third genera- 
tion. The sneutrino mass insertions, which matter for the 
rare decays only to the order we are working, m M = 0, 
are given as (i ^ j) 



8 VL 



1 

r 2 



[cL&USiUsj + x 2 {V VL X L Vl ) tf ) , (24) 



receiving two independent competing contributions. 
Note that the second term in Eq. (24) equals U^S^Uij. 
Since its sign/phase is not fixed by the FN symmetry, 



there can be cancellations in S^f . 



B. A Flavor Model and Phenomenology 

Following Ref. [14] we entertain a U{\) x U(l) FN sym- 
metry with a single symmetry breaking parameter A of 
the order 0.1 — 0.2. Specifically, we use the following 
horizontal charges 



2^(2,0), L 2 (0,2),L 3 (0,2), 
£i(2,l), £ 2 (2,-l),£ 3 (0,-l), 
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M A [GeV] 

FIG. 11: The muon EDM in the flavor model Eq. (25) with 
different fi terms and stau LR-mixings. Fad lines are con- 
sistent with the B(t — > /j/y) bound Eq. (2). In all curves 
tan/3 = 3 and Afj = Mi. The lower two, overlapping curves 
use St LE = 1(T 3 . 



which result in a large 2-3 mixing. The lepton mixing 
angles are obtained as 



A 8 ,(V L Ji 3 ^A 8 ,(F i J 2 



1, 



(Vk)i2 ~ A 2 , (Ve)i3 ~ A 4 , (Ve) 2 3 ~ A 2 , 
<y VL )n ~ f, (26) 

leading to [Zy ~ C(f)- The current level of suppression 
of the 1-3 lepton mixing is considered accidental. 
The slepton soft terms have the flavor structure 

(Xl)i2 ~ A 4 , (Xl)i3 ~ A 4 , (Xl) 2 3 ~ 1, 

(X R ) 12 ~ A 2 , (X R )i3 ~ A 4 , (X R ) 23 ~ A 2 , (27) 

with the diagonal ones {Xl r)h ~ 1. Inserting the fac- 
tors Xl >r and the mixing angles into Eq. (22) the flavor 
changing mass insertions are obtained as 



8 LL " 

°12 


x 2 ,4 X E 

- -A , 6 12 


~^A 2 
r\ 


i °12 




+ x 2 


°13 


x 2 ,4 f.E 


~^A 4 
ri 


i °13 




+ x 2 


°23 


o ^ S E ~ 
ra' 23 


^A 2 , 


°23 


+ 


x 2 ). 



(25) 



The relevant coupling for the muon EDM in this model 
is hence given as 5 2 ^S 23 ~ xix 2 /(r 1 r 2 )\ 2 < 0.04. 

We illustrate the correlation between the rare decays 
and the muon EDM for near maximal mass insertions 
^23 — *0-8 an d S 23 = 0.05 and various higgsino param- 
eters and stau LR-mixings in Fig. 11. Fad lines are in 
agreement with the current B(t — > /17) bound. The two 
upper curves have a negatively valued stau LR-mixing in- 
duced by the /i term. The lower two, overlapping curves 
use S 33 L = 10~ 3 . Here, again, a larger /1 term (dashed 
curve) gives heavier higgsinos and allows for a larger d M 
by relaxing the constraint from the t — > /17 decays. 

The constraint from the FCNC can be avoided in the 
hybrid models by suppressing flavor violation by increas- 
ing the separation between the Planck and the messenger 
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scale. The exact value of the upper bound on x^/r, in the 
model Eq. (25) depends on the susy spectrum and com- 
position and is rather weak due to the strong suppression 
of 1-2 mixing, which passes the fi — > ej constraint. For 
example, for the parameters along the fad lines in Fig. 11 
holds Xi/n < 0(1). 

We obtain in the flavor model Eq. (25) for the maximal 
value of the muon EDM 

^ < 5 • 10- 24 cm, (29) 

e 

in agreement with existing data, specifically the FCNC 
constraints Eq. (2). The other FN models of Ref. [14] 
have smaller values of S^S^a, hence yield smaller val- 
ues of dfj,, either due to a smaller 2-3 mixing, or by a 
stronger constraint on Xi/fi effective for less suppressed 
1-2 mixing. 

We conclude that an observation of the muon EDM 
at the level of order 10~( 23 ~ 24 ) ecm would imply the fol- 
lowing in the context of gauge-gravity models with a FN 
symmetry: 

• The flavor structure needs to be very specific, i.e. 
not every FN symmetry correctly reproducing the 
SM lepton masses and mixings works. 

• CP is broken together with flavor, and the CP 
phases are unsuppressed. 

• The susy spectrum is not too heavy, of the order of 
a few 100 GeV, and y < 1. 

• The messenger scale needs to be very high, not far 
from itim ~ a/(47r)mpi. 

Note that analogous studies in the quark sector have 
found that FCNC data constrain the messenger scale in 
typical FN models to be about three orders of magnitude 
below m pi [15]. 

C. Including ^4-terms 

The gauge mediated contribution to the trilinear 
A-teims is of higher order and usually set to zero 
Ae{™>m) = 0. This has also been employed implicitly in 
the previous section. While the MSSM RG equations do 
induce finite A-terms at the weak scale, however, such 
terms do not introduce CP or flavor violation beyond 
the Yukawa matrices, and do not matter for flavor phe- 
nomenology. 

However, in our hybrid set-up, there can be a finite 
contribution from gravity to the ^4-terms [22]: 

A E (m M ) ~ Vrh^:Y e , (30) 

where we used that the gravity-induced A-terms receive 
the same parametric suppressions from the flavor sym- 
metry breaking as the corresponding Yukawa matrices. 



Here, fh 2 ,x are of the order of rh E and ?7i 2 , x\ and x-i, 
respectively. 

We then evolve according to the MSSM RG equations. 
The result can be recast in the following approximate 
form 

A E {mz) -Y e (a E + Vrh^ + b E Y^Y e ). (31) 

The dominant contribution to as and b E stems from the 
gaugino masses, which drives them to magnitudes of clcc- 
troweak size for large enough m M . Note that a E < and 
b E > 0. Unless tan/3 is very large, the double Yukawa- 
suppression makes the b E term negligible. As indicated 
in Eq. (31), the gravity contribution evolves with a coef- 
ficient of order one. 

In the basis with lepton mass eigenstates and diagonal 
neutralino couplings we obtain: 

A E = V Ll A e V e 

~ D e (a E + Vrh?x + b E D 2 ). (32) 

Under the RG evolution the A-terms mix onto the soft 
masses squared and induce also corrections to the soft 
masses in Eq. (22). Due to the Yukawa texture, this 
results effectively into a correction of the order x of the 
cl,jj coefficients, and will be absorbed therein for the 
purpose of this work. 

The strongest bound on the hybrid model including A- 
terms Eq. (30) is from the electron EDM, stemming from 
E and its imaginary part. An approximate expression 
for d e reads as, see, Eq. (8) for the corresponding formula 
for d M , 

<, M ^y^ ,33, 

where we assumed y = 1 in the second line. Here, CP 
violation is induced by the gravity-mediated ^4-terms, 
Eq. (30), which are fixed by the FN symmetry up to or- 
der one, in general complex numbers, only. One obtains 
from Eqs. (1) and (33): 

x/r < 5 ■ 10~ 3 . (34) 

This bound implies a separation between the Planck scale 
and the messenger scale of about up to four orders of 
magnitude, see Fig. 12. Information on the hadronic sec- 
tor, that is, the neutron EDM in a similar set-up with 
Yukawa-like A-terms and a viable quark FN symmetry 
gives somewhat milder constraints [22]. 

The bound Eq. (34) excludes a muon EDM from flavor 
since the mass insertions are subsequently suppressed as 

^3 S §3 ~ A fc (a;/r) 2 < AMO" 5 < lO" 5 , (35) 

where k denotes the non-negative integer power from the 
FN charges of a given model. 
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FIG. 12: The ratio x/r in minimal gauge mediation Eq. (18) 
as a function of the messenger scale. The curves correspond 
to TVju = 1, 3 and 7 (from top to bottom). The horizontal line 
denotes the constraint from the electron EDM, Eq. (34). 



So far, we assumed the CP phase <fi of Ae n to be 
fully unsuppressed. The bound Eq. (34) weakens by a 
factor 1/sin 2 <f> if <fi is not maximal. However, for k = 2 
as in the model from the previous section, even with a 
phase suppression as low as sin <f> at the percent level, 
the prediction for is as in the linear mass scaling case, 
Eq. (3). 



V. SUMMARY 

Observing a large muon EDM > 0(1O -24 ) ecm re- 
quires CP violation from flavor violation and favors a 
light mass spectrum. The heavier the spectrum and the 
larger the higgsino admixture, the larger is the requi- 
site tuning to avoid the most relevant flavor constraint, 
B(r-»/i7). 

In the context of hybrid gauge-gravity models a mea- 
surement of such large values of the muon EDM would 
point towards a specific Planck scale flavor structure and 
a very large messenger scale not far from ~ a/(47r)mpi. 
In the presence of Yukawa-textured ^4-terms at the level 
of Eq. (30) with CP phases, however, the muon EDM 
is limited by the electron EDM data to follow the lin- 
ear mass scaling constraint Eq. (3), which is below the 
current experimental sensitivity. 

We conclude that the muon EDM is informative on 
various aspects of the underlying fundamental physics. 



Acknowledgments 

KH and TR gratefully acknowledge the support from 
the Academy of Finland (Project No. 115032). The work 
of JL is supported in part by the Foundation for Funda- 
mental Research of Matter (FOM) and the Bundesmin- 
isterium fur Bildung und Forschung (BMBF). 



Appendix A: The relevant MSSM parameters 

The notation follows closely Ref. [10, 20]. 

1. Mixing matrices 

The slepton mass-squared matrices Eqs. (4) and (6) 
are diagonalized by the matrices R 1 '": 



R^Mf^R 1 /*)- 1 = diag(Mj-) 



(Al) 



The neutralino and chargino mass matrices are as usual 

/ Mi 








M 2 



Mv 



9' vi g'vi \ 
gvi gv2 

9' vi gvi 

2^2 _gv2 



M 2 gv 2 
gvi n 



(A2) 



where Mi and M 2 are the U(l) and SU(2) gaugino soft 
masses, respectively, and g' — e/ cosdy/ — y/3/5gi, 
g — e/ sin 9w = 52 ■ The neutralino and chargino mass 
matrices are diagonalized by matrices N and U, V as fol- 
lows: 



WM^V- 1 = diag(M x ±), 
A^*M x oiV _1 = diag(M x o). 

2. Leptonic observables 



(A3) 



The relevant Lagrangian involving sleptons, neutrali- 
nos and charginos is as follows: 



where 



,R = 

-ijk - 

n L = 

hjk = 

i R = 

Hjk — 



C = li(nf jk L + nf jk R)x% 

+k{cf jk L + c? jk R)xJh + h-c, 

-gV sl R%, 

-g'V2N*iRl* i+3 - Vl N* 3 Rt 
(g'Nji + gN j2 ) ifa/yfi - y*N j3 R{* i+3 



(A4) 



(A5) 



The susy contributions to the radiative lepton decay 
width, the lepton EDM, and the magnetic moment can 
be written at one-loop as [10] 



16 

4i Im( " 4 + a " ) ' 
-Re(a l L i +a^), 



d, 
Aa 4 



(A6) 
(A7) 

(A8) 
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where terms of order 0(mi i /mi ) have been dropped. where 



The coefficients a i •' are given as 



4 6 



16 



' m 2 ^ ^ \ \ ikr J fer ^n2 ^ 'hkr^jk 




r 9 



^23 

EE 



/ „L L* m lj , B ff* W M 
WfcrWfcr „9 ^ikr^jkr „9 




_ , . 2 + 3a; - 6a; 2 + a; 3 + 6adoga; , . , „. 
= • ( Al «) 



6(1 -x) 4 
1 - a; 2 + 2a;loga; 
(1 - a;) 3 ' 
1 - 4a: + 3a: 2 - 2a; 2 loga; 

(1 - a:) 3 ' 



(A12) 
(A13) 



The functions G1.3 relevant for the calculation in 
the mass insertion approximation are obtained from 
Eqs. (A10) and (A12) as G 3 (l/x) = (x 2 /2) ■ dF 1 /dx and 
G 1 (l/x)^(x 2 /2)-dF 3 /dx. 
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